EPJ manuscript No. 

(will be inserted by the editor) 



Magnetization of ballistic quantum dots 
induced by a linear-polarized microwave field 

A. D. Chepelianskii 1 and D. L. Shepelyansky 2 

1 Ecole Normale Superieure, 45, rue d'Ulm, 75231 Paris Cedex 05, France 

2 Laboratoire de Physique Theorique, UMR 5152 du CNRS, Univ. P. Sabatier, 31062 Toulouse Cedex 4, France 

October 25, 2006 

Abstract. On a basis of extensive analytical and numerical studies we show that a linear-polarized mi- 
crowave field creates a stationary magnetization in mesoscopic ballistic quantum dots with two-dimensional 
electron gas being at a thermal equilibrium. The magnetization is proportional to a number of electrons 
in a dot and to a microwave power. Microwave fields of moderate strength create in a one dot of few 
micron size a magnetization which is by few orders of magnitude larger than a magnetization produced 
by persistent currents. The effect is weakly dependent on temperature and can be observed with existing 
experimental techniques. The parallels between this effect and ratchets in asymmetric nanostructures are 
also discussed. 

PACS. 75.75.+a Magnetic properties of nanostructures - 73.63.Kv Quantum dots - 78.70.Gq Microwave 
and radio-frequency interactions 



1 Introduction 

Since 1990, when magnetization of an ensemble of 10 
mesoscopic rings had been detected experimentally pQ, 
the problem of magnetization of quantum dots of two- 
dimensional electron gas (2DEG) or persistent currents 
[ attracted a great deal of attention (see e.g. and Refs. 
therein). It is well known that a magnetic field gives no 
magnetization in a classical system at thermal equilibrium 
(see e.g. [1]). Thus, the persistent currents have a quantum 
origin and are relatively weak corresponding to values of 
one electron current of typical strength 3 x 10~ 3 evF/ L p 
where vf is the Fermi velocity and L p is a perimeter 
of quantum dot. Due to that skillful experimental efforts 
had been required to detect persistent currents of strength 
Iq = evp/Lp in an isolated ballistic ring [5]. 

The effects of ac-ficld on magnetization of mesoscopic 
dots have been discussed in [Hid- It was shown that ac- 
driving induces persistent currents which strength oscil- 
lates with magnetic flux. The amplitude of the current is 
proportional to the intensity of microwave field but still 
its amplitude is small compared to one electron current 
evp/Lp. Such currents have pure quantum origin and are 
essentially given only by one electron on a quantum level 
near the Fermi level. Experimental investigations of mag- 
netization induced by ac-driving have been reported in [HI 
|5j . The amplitude of induced currents was in a qualitative 
agreement with the theoretical predictions 017]. However, 
the strength of currents in one ring was rather weak (less 
than nA) and it was necessary to use an ensemble of 10 5 
rings to detect induced magnetization. 



In this work we show that a liner-polarized microwave 
field generates strong orbital currents and magnetization 
in ballistic quantum dots. The magnetization (average non- 
zero momentum) of a dot is proportional to the number of 
electrons inside the dot and to the intensity of microwave 
field. The sign of magnetization depends on orientation of 
polarization in respect to symmetry axis of the dot. It is 
assumed that without microwave field the dot is in a ther- 
mal equilibrium characterized by the Fermi-Dirac distri- 
bution with a temperature T. A microwave field drives the 
system to a new stationary state with non-zero stationary 
magnetization to which contribute all electrons inside the 
dot. The steady state appears as a result of equilibrium 
between energy growth induced by a microwave field and 
relaxation which drives the system to thermal equilibrium. 
In contrast to the persistent currents discussed in P]|^][7| 
the dynamical magnetization effect discussed here has es- 
sentially classical origin and hence it gives much stronger 
currents. However, it disappears in the presence of disor- 
der when the mean free path becomes smaller than the 
size of the dot. Thus the electron dynamics should be bal- 
listic inside the dot which should have a stretched form 
since the effect is absent inside a ring and is weak inside 
a square shaped dot. Also the effect is most strong when 
the microwave frequency is comparable with a frequency of 
electron oscillations inside the dot. The above conditions 
were not fulfilled in the experiments [BE] and therefor the 
dynamical magnetization had not be seen there. 

The physical origin of dynamical magnetization can 
be seen already from a simple model of two decoupled 
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dissipative oscillators for which a monochromatic driving 
leads to a certain degree of synchronization with the driv- 
ing phase |10|. This phenomenological approach was pro- 
posed by Magarill and Chaplik who gave first esti- 
mates for photoinduced magnetism in ballistic nanostruc- 
tures. Here we develop a more rigorous approach based 
on the density matrix and semiclassical calculations of dy- 
namical magnetization. We also extend our analysis to a 
generic case of enharmonic potential inside the dot that 
leads to qualitative changes in the magnetization depen- 
dence on microwave frequency. We also trace certain par- 
allels between dynamical magnetization, directed trans- 
port (ratchets) induced by microwave fields in asymmet- 
ric nanostructures [T^lll3lll4llTS] and the Landau damping 

mm. 

The paper has the following structure: in Section II 
we consider the case of a quantum dot with harmonic po- 
tential, enharmonic potential is analyzed in Section III, 
a quantum dot in a form of Bunimovich stadium |18j is 
considered in Section IV, discussions and conclusions are 
given in Section V. 



2 Quantum dots with a harmonic potential 

Electron dynamics inside a two-dimensional (2D) dot is 
described by a Hamiltonian 

H = (p| + Py) /2m + U (x, y) -xf x cos cut -yf y cos wt (1) 

where m is electron mass and p XtV and x, y are conju- 
gated momentum and coordinate. An external force f xy 
is created by a linear-polarized microwave field with fre- 
quency uj. The polarization angle 9 and the force ampli- 
tude / are defined by relations f x = / cos 6*, f y — fsinO. 
In this Section we consider the case of a harmonic poten- 
tial U(x,y) = m(uj 2 x 2 +ui 2 y 2 ) /2 where oscillation 
frequencies in x, y directions (generally non equal). 

Let us consider first a phenomenological case when an 
electron experiences an additional friction force F = — 7P 
where 7 is a relaxation rate (this approach had been con- 
sidered in and we give it here only for completeness). 
The dynamical equations of motion in this case are linear 
and can be solved exactly that gives at t ^> 1/7: 



x{t) = SR-I: 



t -J*lH— — 3? Y(t), 



■u) 2 -\-i^u) 



(2) 



where 3? marks the real part. Then, the electron velocities 
are 

x (t) = ft /m _ 



y (t) = ft te a , u /y( m - 



= »y(t). 



(3) 



This gives the average momentum 

L = m (x(t)v y (t) - y(t)v x (t)) 
-iwfxfy/m 



n 



ijuj) 



(4) 



In the limit of small 7 -C 0J x , y Eq.(0} gives for the off- 
resonance case 



"off 



7 w 2 (o)2 _ uj 2 )f x f y /m 



(uj 2 -^) 2 ^ 2 -^) 2 - 



while at the resonance lu = uj x 



fxfy/m 
li^l - u 2 ) ' 



(5) 



(6) 



From a physical viewpoint an average momentum appears 
due to a phase shift between oscillator phases induced 
by dissipation and an orbit takes an elliptic form with 
rotation in one direction. In some sense, due to dissipation 
the two oscillators become synchronized by external force 

. As usual 0] , an average orbital momentum L for one 
electron gives a total magnetic moment M = NLe/2mc 
where N is a number of electrons in a quantum dot. A 
pictorial view of spectral dependence of M on to at various 
ratios to x /io y is given in 

To extend the phenomenological approach described 
above we should take into account that the electrons inside 
the dot are described by a thermal distribution and the ef- 
fects of microwave field should be considered in the frame 
of the Kubo formalism for the density matrix (see e.g. [3]). 
For analysis it is convenient to use creation, annihilation 

operators defined by usual relations x = \J 2 mu> (° + 



£-(6 + b+) and p x = -1 



2 m lu 



ihujj- 



(a - a+), py 



The unperturbed Hamiltonian in ab- 
sence of microwave field takes the form Hq = htu x (a + a + 
1/2) + Tiujy(b + b + 1/2). Then the orbital momentum is 



L = 



if l 



2yjLO X Uly 



[(id x — ui y )(ab — a + b + ) + 
(ui x +Ld y )(ab — a b)]. 



(7) 



The only non zero matrix elements of L are 
< n x ,n y \L\n x + S x , n y + 5 y > = 

: (5 X U) X — SyLOy) X 

2^LU X 0Jy 

[(n x + (1 + S x )/2)(n v + (1 + S v )/2)}^ 2 , (8) 

where 8 XtV = ±1 and n xy are oscillator level numbers. 
The perturbation induced by a microwave field V(t) — 
{—fxX — f y y)cosujt — Vcoscot should be also expressed 
via operators a, a + , b, b + . 

In the Kubo formalism the evolution of the density 
matrix p(t) is described by the equation 



ih 



HA 



[H + V{t),p]-itvy{p-po) 



(9) 



where po is the equilibrium density matrix: 
Pa = J2n x ,n y >o Pr^,n y \n x ,n y >< n x ,n y \. Using perturba- 
tion theory p(t) can be expanded in powers of the external 
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potential amplitude V(t) : p(t) = p + pi{t) + fait) + < n x + 5 X , n y \pi(t)\n x ,n y >>t 



In the first order we have 



<< n x + S x , n y + 6 y \pi(t)\n x ,n y + 8y > x 



< a|pi(i)|/3 >= Z p ^ a ^ < a \ V \P > c ^t/2 | < n x,n y +S y \V(t)\n x ,n y » t = 

e - e a - fiuj + ijh - [(n lJ + (l + ,5,)/2)K + (l + 4)/2)/(64m 2 ^)] 1 / 2 



(p(ep) - p(e a )) < a\V\(3 > e _ iu , x 1 , , 



-e~ tuJt /2 . (10) 

For the harmonic potential we obtain (P™x,n„ — Pn^+s^^y — Pn xi n v +s y + Pn^+s^^y+Sy) 

< n x + S X , n y \pi(t)\n x ,n y >= and 

, w /n a; + (l + 4)/2 <n x + 4,% + 5j / | </3 2 (i) >t K>% >= 

(Pn x ,n H - Pn x +5 x , n y)jx\ X (11) 



%Timuj x fxf y 



.., -9n a ,n v ,S x ,5 v (w) X 

w x 5 x - lu + 27 + -lu x 5 x +lu + ij J ' [(n x + {l + 5 x )/2)(n y + (l + 5 y )/2)]^ 2 , 

where 

< n x ,n y +S y \p 1 (t)\n x ,n y >= y = 

)fv\ — ^x (12) 



Shmujy —u) x 5 x —07 + 27 — Wj^j; + 0^ + 27 

1 1 



-WyOy — + 27 —<jJ y o y + Ld + J y y ' V V < ' I 

Pn x ,n y ~ Pn x +8 x ,n y ~ Pn x ,n y +8 V + Pn^+S^ ,n y +5 y 

The time averaged second order correction to the density —ui x 8 x — ui y 5 y + 27 

matrix is given by 

Thus 

<a\< p 2 [t) > t |/3 >= tp -e a + ifry < n -' n 2/l X K + n v + S v > x 

(13) < ?i K + 5^, n y + 5„| < /3 2 (t) > t \n x , n y >= 

where [...] marks the commutator between two operators. if x f y 

To compute the average momentum L we need to find lQmcu to ^ n *> n y> n *+ s *> n y+ s y 

the terms < n x + 6 x ,n y + S y \ < p2(t) >t \n Xl n y >. Ac- 
cording to (fH5f they are expressed via the matrix elements where 
like 

« n x + S X: n y + S y \V(t)\n x + S X: n y > x 1 1 

<n x +d x ,n y \ Pl (t)\n x ,n y » t = (14) ^ -uj x S x - uj + 27 + -uj x 5 x + uj + 27 + 

).fxfy X 1 _^ 1 x 

[(n. y + (1 + Sy)/2)(n x + (1 + 4)/2)/(64m 2 w x ^)] 1 / 2 x - w + £7 + w + £7 

1 1 \ S X UJ X — SyWy 



and 



-W^Jx - UJ + 27 -CJ X (5 X + W + 27/ ' -^xfe - Wj/^ + 27 

(Prix,n y ~ Prix+Sx^iy ~ Pn^^y+Sy + Pn x +S xt n v +5y) x 



<< n x + 5a,, n y + 5j / |pi(t)|n x ,n J/ + 5 y > x 
< n x ,n y + S y \V(t)\n x ,n y » t = (15) 

{Pn^^iy+Sy — Pn x +5 x ,n y +8 y )fxf y x 

[(n, + (l + ^)/2)(n :c + (l + 5 :c )/2)/(64m 2 a; ;E ^)] 1 / 2 x = : G nx ,„ y ,„ x+5x ,„ H+5H (w) 

LomujxUJy , 

1 1 \ n K ,n y ^u,^^ ,o y 



K + (1 + 0/2) (n„ + (1 + 5„)/2) . 
Therefore, the final result is 

< L >= Tr(L < p 2 (t) >) 



uj x S x — uj + 27 — lu x S x + uj + 27 J _ fxfy _ 



Ybmu) x LUy — ' 

Therefore »x,«„>o 



_ fxfl) j/ \ j^j 

« n x + S x , n y + 5 y \V(t)\ n x + S x , n y > x 16muj x u> y ' 
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where 



1 



1 



uj x 5 x - uj + i-f -uj x 5 x + uj + i-f 



1 



1 



-IjJySy — UJ + «7 -WySy + UJ + 1^ 

S X S V . 



(22) 



S X <^X — SyOJy 



-UJ X 5 X - UJySy + 17 



Here N is the number of electrons in the quantum dot. 
Of course, I{ui) is real and can be presented by another 
equivalent expression: 



R(uj) ' 



(23) 



with 



Q(u) = 8-foj x uj y (LJ x - uj 2 ){5j b + + 9 7 2 w 4 + 



2uj 6 + (6 7 4 + 15 7 2 w 2 - uj 4 )uj 2 y + 

(j 2 + 3uj 2 )uj y +uj x ( 7 2 + 3c 2 +^) + 

w 2 (6 7 4 + 15 7 2 cj 2 - uj 4 + (7 7 2 - 8uj 2 )ujI + uj 4 )) 

and 

R(uj) = (( 7 2 + uj 2 ) 2 + 2( 7 2 - uj 2 )uj 2 x + ujt) x 

((7 2 +c 2 ) 2 + 2(7 2 -c 2 )c 2 +^)x 

(cut + uj y + 7 4 + 2 [(co 2 + c 2 ) 7 2 - uj 2 x uj v ]) . 

We used the relation 



(24) 



(25) 



(n x + (1 + 6 x )/2)(n y + (1 + <5. y )/2) + 

(n, + (1 - 5 x )/2)(n y + (1 - 5 y )/2) - 

(n x + (l-S x )/2)(n y + {l + 5 y )/2)- (26) 

(n x + (1 + 4)/2)K + (1 - 5 y )/2) = S x S y 

to reduce Eqs.JTnjl,® to Eas. (|2T)) - ijg3|) . 

It is important to note that the final result l|22|l for the 
average momentum L is independent of unperturbed ther- 
mal distribution p na: ,n y - The momentum grows linearly 
with the number of electrons in the dot N. For uj ~ 7 ~ 
uj x ~ uj y we have Z/ ~ fxfyN/fauj 3 ) in agreement with 
the phenomenological result However, the exact de- 
pendence H22l) - I(25() obtained here from the Kubo theory 
is different from the phenomenological result (0J obtained 
originally in For example, at uj 3> uj x ~ uj v > 7 
our result gives L ~ jN/(muj 2 (uj 2 — w 2 )) while the phe- 
nomenological result (0) of Jl] gives £ cx 1/cj 6 . 

To obtain the expression for L we used above the quan- 
tum Kubo theory. However, the result 1)22(1 has a purely 
classical form and therefore it is useful to try to obtain it 
from the classical kinetic theory. With this aim let us con- 
sider an arbitrary two dimensional system with a Hamil- 
tonian H = H(q x ,p x ,q y ^p y ,t). Then the kinetic equation 
for the distribution function p(x) JJ] reads 



dp 
Of 



+ {H,p} = -7O-P0) 



(27) 




Fig. 1. (color online) Dependence of the rescaled magneti- 
zation M r — —2M/Mo on the rescaled microwave frequency 
Lojux for a quantum dot with a harmonic potential at uj v /uj x — 
3 (see definitions of Mo and M r in Eq 13311 ). Here L x = vp/u x 
and vp is the velocity Fermi at the Fermi energy Ep. The 
rescaled relaxation rate is 7/0^ = 0.1; 0.25; 0.5 (red, blue, green 
curves/circles from top to bottom at uj/u) x = 1). Circles show 
numerical data obtained from Eq. I30H by integration of classi- 
cal dynamics and Monte Carlo averaging over 10 4 trajectories 
from the Fermi-Dirac equilibrium distribution at zero temper- 
ature; full curves show the theoretical result given by Eas.iSTt- 



where po is the equilibrium thermal distribution, {iJ, /} 
are the Poisson brackets and for simplicity of notations 
x = (x,p x ,y,p y ). After a change of variables (i',x') = 
(t, To,* x ) this equation is reduced to 



dp(x',t') 
dt' 



- 7 (p(x',i')-Po(lV,ox')) (28) 



where Tf t t notes the time evolution operator from time t 
to time t' given by the dynamics of the Hamiltonian H. 
This equation can be solved explicitly that leads to the 
time averaged distribution function : 

r° 

< p(x,t) > t = 7 / e T* < /5o(T t+t , )t x) > t dP (29) 



where Tt+t',t is the time evolution operator in the phase 
space. Then the average of any quantity Q(x) (for example 
Q = xp y — p x y) can now be expressed as : 

< Q(x) >= J Q(x) < p(x,t) > t dx = 

7 < [ ( Q(x)p (T t+t ,, t x)dx)dP > t = 

J — OO J 

7 I dx Po (x) [ dPe^' < Q(T w x) > t = 

J J — OO 

7 J dxpo{x)J dPe~^' < Q(T t+t ,, t x) > t . (30) 
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Here it is used that the transformation x' = T t+t ' jt x is 
area-preserving in the phase space. For the case of two 
oscillators with frequencies cj x ,lu v the time evolution can 
be find explicitly so that for dynamics in x we have 



'x(t + t') 



cosuj x t' sintu x t' 
— sin.cj.7-t' coscj T t' 



— h-l^L cos uit 

til* — M* 5 



■ sin tot 



(31) 



with a similar equation for y,p y . After averaging over t we 
obtain 



fx 



fy/m 



<x(t + t') Py (t + t') >^2{^-^)(^-u 
[ sm(uj y t') cos(u! x t') + smut' cos(uj x t') + 

iO 

— sm{uj x t ) cos(ui y t ) smfwjt ) coswt J 



(32) 



with a similar expression for < y(t + t')p x (t + t') > t . After 
substitution of lj3*2)l in l|3TJ|) the integration over t' gives ex- 
actly the expression with I(u>) given by Eqs. (J23Il- (E5J) - 
The integration can be done analytically or with help of 
Mathematica package. This result shows that the average 
momentum L can be exactly obtained from the classical 
formula (|3T)jl . 

A comparison of results of numerical simulations of 
classical dynamics with Monte Carlo averaging over large 
number of trajectories from an equilibrium distribution 
is shown in Fijfl] The numerical data clearly confirm the 
validity of the theoretical expressions given by Eqs.JS^J- 

It is interesting to note that if instead of Eq. @ for the 
density matrix one would assume an adiabatic switching 
of microwave field with a rate 7 then the average induced 
momentum L would be given by the classical relation 0} 
for a classical oscillator. Indeed, such a procedure simply 
induces an imaginary shift in driving frequency. Such type 
of switching had been assumed in and may be con- 
sidered to give a qualitatively correct result even if a rig- 
orous description is given by Eq.© with the final answer 
in the form of Eas. 1)21125(1 being quantitatively different 
from Eq. © . 

For comparison with the physical values of magnetic 
moment M in real quantum dots it is convenient to use 
rescaled momentum M r . To do this rescaling we note that 
the magnetic moment is expressed via the orbital mo- 
mentum as M — eL/(2mc) (in SGS units). Due to the 
relation l|21|l it is convenient to choose a unit of orbital 
momentum induced by a microwave field for one electron 
as I/o = mv FL x (f x f y L 2 )/E F where Ep — mv F /2 is the 
Fermi energy in a dot. Then the unit of magnetic momen- 
tum is Mo = eiVLo /mc where N is the number of electrons 
in a dot. This implies that the physical magnetic moment 
M can be expressed via our rescaled value according to 




Fig. 2. (color online) Dependence of rescaled momentum M r 
given by Eg, 1331 on rescaled frequency ujL x /vf for the case 
of quartic oscillator with the Hamiltonian 1341 at K — and 
L y /L x = l/\/3 {r x /r y — 1/3), the monochromatic field has 
fx = f y and f x L x /EF = 0.25??. Curves correspond to differ- 
ent values of the relaxation rate *yL x /vF — 0.03 (red), 0.06 
(blue), 0.14 (green) (from bottom to top at u>L x /vf = 1.2). 
The curves are obtained from Eq.J^OJ by integration of classi- 
cal dynamics and Monte Carlo averaging over 10 4 trajectories 
from the Fermi-Dirac distribution at zero temperature. 



the relations 



M = -M r M /2, M = eNLo/mc, 
M = eNv F L x (f x f y L 2 x )/(cE 2 F ). 



(33) 



where the oscillation frequency in x is lj x = Vf / L x . It is 
important to stress that the total magnetization of the dot 
is proportional to the number of electrons in a dot with 
fixed E F . 



3 Dots with an enharmonic potential 

It is very important to extent the methods developed 
above to a generic case of enharmonic potential inside a 
dot. With this aim we consider the case of 2D quartic 
nonlinear oscillator described by the Hamiltonian 

H=(pl+ P 2 y )/2m + (r 2 x x 4 + r 2 y y 4 + 2Kx 2 y 2 )/4. (34) 

For K = we have two decoupled quartic oscillators. Due 
to nonlinearity the frequencies of oscillations scale with 
energy as oj X:V ps 1.2(r 2 Ep/m 2 ) 1 ^ (see e.g. |HJ). In this 



case L x = {AE F /rl) x / A , L y 



(AEp/r 2 ) 1 / 4 and we choose 
1/V3 for our stud- 



the fixed ratio L y /L x — \Jr x jr y 
ies. The rescaled momentum and magnetization are again 
given by Ea. (|3*3l . 

The value of averaged orbital momentum is obtained 
by numerical integration of hamiltonian equations of mo- 
tion in presence of a linear-polarized monochromatic field 
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Fig. 3. (color online) Same as in Fig. 2 but for ^L x /vf = 0.14 
(green), 0.3 (violet), 0.45 (yellow), 0.55 (black) (from top to 
bottom at u)L x /vf = 1.2). 

(see (JsHJl). The effects of relaxation to stationary state 
with a rate 7 are taken into account via relation (|3L)[) . The 
average momentum < L > is obtained by Monte Carlo av- 
eraging over 10 4 trajectories from the Fermi-Dirac distri- 
bution at zero temperature. The integration time is about 
10 4 oscillation periods. 

The dependence of rescaled magnetization on microwave 
frequency 10 for different relaxation rates is shown in FigsEEl 
In contrast to the case of harmonic potential the depen- 
dence on frequency is characterized by a broad distribu- 
tion with a broad peak centered approximately near os- 
cillation frequency lu x — \.2(r x Ep jm 1 ) 1 ^ w \.2vp/L x . 
Significant magnetization is visible essentially only inside 
the interval uj x « 1.2vf/L x < lu < uj y w 2.\vp / L x . Data 
obtained also show that at small relaxation rates mag- 
netization drops to zero approximately as M r cx 7 (see 
Fig [21). Also M r drops with the increase of 7 at large 7 
(see Fig|3J|. This behaviour is qualitatively similar to the 
case of harmonic potential (see Eqs.©,©)- 

It is important to note that according to data shown 
in Fig^the coupling between x, y-modes, which generally 
leads to a chaotic dynamics does not lead to signif- 
icant modifications of the magnetization spectrum. Only 
at rather large values of K, when the modes are strongly 
deformed, the spectrum starts to be modified. 

To show that the magnetization dependence on fre- 
quency found for the quartic oscillator (fMj) represents a 
typical case we also study magnetization in chaotic bil- 
liards described in the next Section. 



4 Magnetization in chaotic billiards 

To study microwave induced magnetization in billiards we 
choose the Bunimovich stadium billiard ^B] as a typical 
example. The semicircle radius is taken to be R, the total 
size of stadium in x is L x , and in y it is L y = 2R. Usually 



0.3 




Fig. 4. (color online) Same as in Fig. 2 for ^L x /vf = 0.14 
and different coupling strength K between x,y- modes in (13411 : 
KL%/E F = 4 (red), (green), 8 (blue), 16 (black) (from top 
to bottom at uiL x /v F = 1.3). 




I , 1 1 1 1 1 1 

-2-1012 
x 

Fig. 5. (color online) Example of a trajectory inside the Buni- 
movich stadium dot in presence of microwave driving. The 
system parameters are L x /R — 3.5, L y /R = 2, microwave 
polarization angle 6 = 7r/4 (f x = f y ), fR/E F — 0.28, 
loL x /V f = 1.7, T/E F = 0.1, v F r rel /R w 2, v F n/R w 500. 

we use L x = 3.5R (see Fig|5J). Inside the billiard the parti- 
cle is affected only by monochromatic force, the collisions 
with boundaries are elastic. 

To take into account that without monochromatic force 
the particles relax to the Fermi-Dirac equilibrium we used 
the generalized Metropolis approach developed in |14j . 
Namely, after a time interval At the kinetic energy of 
particle E is changed randomly into the interval (E — 
AE, E + AE) with a probability given by the Fermi-Dirac 
distribution at given temperature T. The change takes 
place only in energy while the velocity direction remains 
unchanged (see below). This procedure imposes a conver- 
gence to the Fermi-Dirac equilibrium with the relaxation 
time r re i w At(Ep / AE) 2 (see Jl] for detailed description 
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Fig. 6. (color online) Rescaled momentum M r as a function 
of rescaled microwave frequency u)L x /vf in the Bunimovich 
billiard. All parameters are as in Fig. 5, curves correspond to 
different values of relaxation time r re ( with VFT re i/R = 0.5 
(green), 0.7 (violet), 1.1 (blue), 2 (red), 8 (black) (from top to 
bottom at ujL x /Vf = 2). 




2 3 

coL x / V F 

Fig. 7. (color online) Same as in Fig[(j]but the length of sta- 
dium is increased so that here L x /R — 5. 



of the algorithm). In the numerical simulations we usually 
used vpAt/R « 0.05 and AE/E F = 0.15. At such pa- 
rameters a particle propagates on a sufficiently large dis- 
tance during the relaxation time: vpT re ijR ~ 2. In absence 
of ac— driving the Metropolis algorithm described above 
gives a convergence to the Fermi-Dirac distribution with 
a given temperature T . In presence of microwave force the 
algorithm gives convergence to a certain stationary distri- 
bution which at small force differs slightly from the Fermi- 
Dirac distribution (the dependence of deformed curves on 
energy is similar to those shown in Fig|2Jof Ref. How- 
ever, in contrast to the unperturbed case, the perturbed 
stationary distribution has an average nonzero orbital and 
magnetic momentum M. The dependence of average mo- 
mentum M on temperature T is relatively weak if T <§C Ep 
and therefore, the majority of data are shown for a typ- 
ical value T/Ep — 0.1 (see more detail below). To take 
into account the effect of impurities the velocity direction 
of a particle is changed randomly in the interval (0, 2n) 
after a time r,. Usually we use such Tj value that the 
mean free path is much larger than the size of the billiard 
VFTi/R 500, in this regime the average momentum is 
not affected by Tj (see below) . The average momentum is 
usually computed via one long trajectory which length is 
up to 10 times longer than R; computation via 10 shorter 
trajectories statistically gives the same result. An example 
of typical trajectory snapshot is shown in FigOO It clearly 
shows a chaotic behaviour (the lines inside the billiard are 
curved by a microwave field). 

The numerical data for dependence of average rescaled 
momentum M r (see Eq. ((3SJ) on rescaled microwave fre- 
quency u)L x /vp are shown in Fig. for polarization 9 = 
7r/4. Qualitatively, the dependence is similar to the case 
of nonlinear oscillator discussed in previous Section (see 
Figs. 121 II) . At the same time, there is also a difference in 




I i i i i i i i 

-2-1012 
x 

Fig. 8. (color online) Example of a trajectory inside a dot in 
a form of the Bunimovich stadium with a circular "impurity" 
inside the billiard. The billiard boundary is as in Fig0 the 
circular impurity has the radius r = R/2, its center is located 
at x — —0.75R, y — R/3 counting from the center of the 
stadium. The system parameters are as in FigU3 except that 
the polarization angle 6 — 0. 

the behaviour at small frequency (u)L x /vp < 1.5) where 
the momentum changes sign. At small relaxation times 
T re i the frequency dependence has a sharp peak near 
ujL x /vf « 2. The increase of r re i leads to a global de- 
crease of average magnetization, that is similar to the data 
of Fig. [3J also the peak position shifts to a bit higher u>. 
Let us also note that according to our numerical date the 
rescaled momentum M r is independent of the strength of 
driving force / in the regime when fR/Ep < 0.5. This is 
in the agreement with the relation l|33|l . 

The position of peak is determined by the frequency 
of oscillations along long x-axis of the billiard. Indeed, an 
increase of this size of billiard from L x = 3.5i? (Fig. [SJ 
to L x = 5R (Fig. keeps the shape of resonance curves 
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Fig. 9. (color online) Dependence of rescaled momentum 
M r ~ —2M/Mq on rescaled frequency loL x /vf for the Buni- 
movich stadium with impurity shown in Fig|5] Here /R/Ef = 
0.28, 9 = (bottom black curve) and fR/E F = 0.2, 9 = vr/4 
(top red/gray curve); also in this Figure we use definition 
Mq = eNvpL x (f 2 L x )/(cE F ) which is more suitable for po- 
larization with 9 — 0. Other parameters are T/Ef = 0.1, 
VFTrei/R ~ 2, VFTi/R « 500 (cf. with corresponding case in 
Fig.EJ. 




Fig. 10. (color online) Dependence of momentum M on polar- 
ization angle 9 (M is rescaled to its value at 9 — 7r/4) for the 
Bunimovich stadium (red/gray) and the stadium with circular 
impurity (black). Numerical data are shown by open circles 
obtained for ujL x /vf = 0.78 (red/gray) and uiL x /vf = 1-56 
(black) at fR/E F = 0.28, T/E F = 0.1, v F r re i/R « 2, 
vf^/R « 500. Smooth red/gray curve shows the theoret- 
ical dependence sin 29; black curve shows a numerical fit 
M{9)/M(tt/4) = -0.15 + 1.23sin(26> - 0.32). 



practically unchanged. At the same time the rescaled mag- 
netization drops approximately by a factor 2. This means 
that there is no significant increase of M with increase of 
L x (Mq grows by a factor 2.9). From a physical view point, 
it is rather clear since in the regime with L x 3> L y further 
increase of L x cannot lead to increase of magnetization. 
This means that in Eq. i|33|) the value of Mo gives a correct 
estimate of real magnetization assuming that L x ~ L y , 

The Bunimovich billiard has varies symmetries, namely 
x — » — x, y — > —y. It is interesting to study the magneti- 
zation properties when all of them are absent. With this 
aim we introduced an elastic disk scatterer inside the bil- 
liard as it is shown in Fig. |SJ The dependence of rescaled 
magnetization M r on rescaled frequency loL x /vf is shown 
for this "impurity" billiard in Fig. for two polarization 
angles 9 — 7r/4 and 9 = 0. For 9 = ir/4 the behaviour in 
Fig-EJis rather similar to the case of billiard without impu- 
rity (see Fig. HJ) , taking into account that Mo values differ 
by a factor two that gives smaller M r values in Fig. In 
addition the peak at u}L x /vp — 2 is more broad that can 
be attributed to contribution of orbits with a shorter pe- 
riods colliding with the impurity. The case of polarization 
with 9 = is rather different. Indeed, here the average 
magnetization exists even if in the billiard case M r = 
due to symmetry. Also the sign of the momentum (magne- 
tization) is different comparing to the case of 9 = 7r/4 po- 
larization. For the impurity billiard we find that at 9 = 
the absolute value of magnetization decreases with the in- 
crease of relaxation time T re ; in a way similar to one should 
in Figs. Ed 



The polarization dependence of magnetization is shown 
in more detail in Fig. 1101 for the Bunimovich stadium 
(Fig. EJ) and the impurity billiard (Fig. EJ). In the first 
case we have M(9) oc sin26' as in the case of oscillator 
so that the magnetization averaged over all polarization 
angles is equal to zero. In contrast to that in the second 
case when all symmetries are destroyed the averaging over 
all polarization angles gives nonzero magnetization of the 
dot. In addition to that internal impurity gives a phase 
shift in the polarization dependence. The phase shift is 
due to absence of any symmetry. In such a case we have a 
more general dependence M oc (f x f y + af x + bfy), where 
o, b are some constants. 

In addition, we also checked that if at 9 = the disk 
impurity inside the billiard is replaced by a semidisk of 
the same radius then the dependence on the parameters 
remains essentially the same, as well as the sign of mag- 
netization (the disk is divided by a vertical line onto two 
semidisks and left semidisk is removed). It is interesting 
to note that a negative sign of M r means anti-clockwise 
rotation. This direction of current rotation can be also un- 
derstood in a link with ratchet flow on the semidisk Galton 
board studied in |12U13U14| . The link with the ratchet ef- 
fect becomes especially clear if to consider the semidisk 
ratchet billiard shown in Fig. El Here, due to the ratchet 
effect discussed in |12U13U14| for polarization 9 = elec- 
trons should move to the left in the upper half of the 
billiard and to the right in the bottom half, thus creating 
negative magnetization M r . We think that if only upper 
semidisk is left still the direction of rotation at 8 = 
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Fig. 11. Example of a billiard with ratchet effect. Collisions 
with the stadium and internal boundaries are elastic. 



microwave polarization is due to the particle flow directed 
from right to left in the upper part of the billiard that cor- 
responds to negative sign of magnetization M r in Fig. [!|] 
Finally, let us make few notes about the dependence 
of magnetization M r on the impurity scattering time ti 
and temperature T. For example, for the impurity bil- 
liard of Fig. |H1 at 6 = we find that the variation of M r 
is about 10% when v F TijR is decreased from 500 down 
to 50, and then M r drops approximately linearly with 
v f tJR (e.g. M r (v F Ti/R = lO)/M r (v F n/R = 50) = 0.47, 
M r (v F n/R = 5)/M r (v F n/R = 50) = 0.2, M r (v F n/R = 
2.5) /M r (v F Ti/R = 50) = 0.08. The decrease of magne- 
tization with decrease of 7j is rather natural since the 
ballistic propagation of a particle between boundaries dis- 
appears as soon as the scattering mean free path v F Ti 
becomes smaller than the billiard size L x . Similar effect 
has been seen also for the ratchet flow on the semidisk 
Galton board (see Fig. El in Ref.[T3|). As far as for depen- 
dence on temperature our data show that it was relatively 
weak in the regime T/E F <C 1. For example, for the set 
of parameters of Fig. EH at 8 = the value of M r is in- 
creased only by 14% when T is decreased by a factor 4 
(from T/E F = 0.1 to T/E F = 0.025) and it is decrease by 
40% when T increased from T/E F = 0.1 to T/E F = 0.2. 
This dependence on T shows that the magnetization of 
the dot remains finite even at zero temperature. This be- 
haviour also has close similarity with the temperature de- 
pendence for the ratchet flow discussed in |14U15j . Such 
ballistic dots may find possible applications for detection 
of high frequency microwave radiation at room temper- 
atures. Indeed, the ratchet effect in asymmetric nanos- 
tructures, which has certain links with the magnetization 
discussed here, has been observed at 50GHz at room tem- 
perature [2"T]. 



5 Discussion 

According to the obtained numerical results (see e.g Fig[5J) 
and Ea. (|33[) the rescaled magnetization of a dot can be 
as large as M r /2 ^0.1 when the relaxation time T re i is 



comparable with a typical time of electron oscillations in 
a dot. Thus, the magnetization of a dot is M = O.lMo. 
According to Eq. l|33|l . at fixed electron density the mag- 
netization is proportional to the number of electrons N 
inside the dot. Due to this the magnetization induced by 
a microwave field can be much larger than the magnetiza- 
tion induced by persistent currents discussed in 0|2 3 , 6 
|7]. For example, for 2DES in AlGaAs/GaAs the effective 
electron mass m = 0.067m e and at electron density n e w 
2 x 10 n cm~ 2 we have E F « lOOif, v F /c w 1.4 x 10~ 3 . 
Hence, according to Eq. (|33|l for a microwave field of f/e = 
IV/ cm acting on an electron in a dot of size L x = 1/xm 
we obtain N k 2x 10 3 , JR/E F k, 0.01 and M /fi B ~ 
2m e v F L x N(fL x /E F ) 2 /h » 5 x W 3 N(fL x /E F ) 2 « 10 3 , 
where [Ib = eh/(2m e c) is the Bohr magneton. Nowadays 
technology allows to produce samples with very high mo- 
bility so that the mean free path can have values as large as 
few tens of microns at 4K . At n e = const we have the scal- 
ing Mo cx L^f 2 and for an increased dot size L x — lOfim 
and field f/e = iV/cm the magnetization of one dot is 
M ss 0.1M i=i IOVb being comparable with the total 
magnetization of 10 7 rings in p0|. Therefore, this one dot 
magnetization induced by a microwave field can be ob- 
served with nowadays experimental possibilities [8TI31I27)] . 
We note that the magnetization is only weakly dependent 
on temperature but the mean free pass at given tempera- 
ture should be larger than the dot size. 

This ballistic magnetization should also exist at very 
high frequency driving, e.g. THz or optical frequency, which 
is much larger than the oscillation frequency. In this regime 
the magnetization drops with the driving frequency M oc 
1/uj 2 (see Eqs. I|21 fl -H25 fl ) but this drop may be compen- 
sated by using strong driving fields. 

From the theoretical view point many questions re- 
main open and further studies are needed to answer them. 
Thus, an analytical theory is needed to compute the mag- 
netization in dots with enharmonic potential or billiards. 
On a first glance, as a first approximation one could take 
analytical formulas for harmonic dot Eqs. I|21 |l -I|25 |l and 
average this result over frequencies variation in an en- 
harmonic dot. However, in the limit of small relaxation 
rate 7 (or large r re i) such an averaging gives finite mag- 
netization independent of 7. Indeed, in analogy with the 
Landau damping |lrjU17| such integration gives effective 
dissipation independent of initial 7 — > 0. Appearance of 
such magnetization independent of 7 would be also in a 
qualitative agreement with the results obtained for ratchet 
transport on the semidisk Galton board |13U14| and in 
the asymmetric scatterer model studied in ^2J- Indeed, 
according to these studies the velocity of ratchet is inde- 
pendent of the relaxation rate in energy (relaxation over 
momentum direction is reached due to dynamical chaos). 
In fact this indicates certain similarities with the Landau 
damping where the final relaxation rate is independent of 
the initial one. Also, the results obtained in [141115) show 
that the ratchet velocity can be obtained as a result of 
scattering on one asymmetric scatterer. Therefore, it is 
rather tentative to use the semidisk billiard of Fig. II II and 
to say that the magnetization in it appears as the result 
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of ratchet flow: for polarization 9 — the ratchet flow 
goes on the left in the upper part of the billiard and on 
the right in the bottom part. The velocity of such ratchet 
flow v f in this 2DES is given by the relations found in [141 
ITS] , namely Vf/vp ~ {fR/Ep) 2 . This flow gives a magne- 
tization of billiard dot induced by a microwave field which 
is of the same order as Mq in the relation However, 
in this approach the rescaled momentum M r is simply 
some constant independent of the relaxation rate in en- 
ergy, as it is the case for the ratchet transport on infi- 
nite semidisk lattice. This result is in the contradiction 
with our numerical data for magnetization dependence on 
the relaxation rate in energy which gives approximately 
M r oc jLx/vf ~ L x /(vpT re i) (see Figs l2l6l7l) . A possible 
origin of this difference can be attributed to the fact that 
the ratchet flow is considered on an infinite lattice while 
the magnetization takes place in a confined system and 
the relaxation properties arc different in these two cases. 

Also we should note that here we used approximation 
of noninteracting electrons and neglected all collective ef- 
fects. In principle, it is well known that microwave radi- 
ation can excite plasmons in 2DES (see e.g, [22]). These 
excitations can be viewed as some oscillatory modes with 
different frequencies uj x ,uj y and therefore in analogy with 
Eqs. (|21[l - 125[) it is natural to expect that a liner-polarized 
radiation can also create rotating plasmons with finite 
magnetization induced by this rotation. In addition, the 
effects of screening should be also taken into account. 
All these notes show that further theoretical studies are 
needed for a better understanding of radiation induced 
magnetization in 2DES dots. 

We thank Kvon Ze Don for useful discussions and for 
pointing to us Ref. 11 at the final stage of this work. This 
work was supported in part by the ANR PNANO project 
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